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Abstract 

In this paper we introduce a class of constraint logic programs such that their termination can be 
proved by using affine level mappings. We show that membership to this class is decidable in poly- 
nomial time. 
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1 Introduction 

Termination is well-known to be one of the crucial properties of software verification. 
Logic programming, and more generally constraint logic programming (CLP), with their 
strong theoretical basis lend themselves easily to termination analysis as witnessed by a 
very intensive research in the area. 

In this paper, which is a revised version of ( Serebr enik and Mesnard 20 04 ). we study 
decidability of termination for CLP(C) programs for a given constraint domain C. In 
general, decidability depends on the constraint domain C. On the one hand, Devienne 
et al. (1993 ) have established undecidability of termination for one-rule binary CLP(H) 
programs, where H is the domain of Herbrand terms. On the other hand, Datalog, i.e., 
logic programming with no function symbols, provides an example of a constraint pro- 
gramming language such that termination is decidable. We note that the decidability of 
the related problem of boundedness for Datalog queries has been studied, for instance, 
in (lAfrati et al. 2 005 ; Marcinkows ki 1996j l. For constraint domains with the undecidable 
termination property, we are interested in subclasses of programs such that termination is 
decidable for these subclasses. A trivial example is the subclass of non-recursive programs. 

We organise the paper as follows. After the preliminary remarks of Section [2] in Sec- 
tion|3]we present our main result. Section previews related results before our conclusion. 
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2 Preliminaries 

For CLP-related definitions, we follow (Jaff ar et al. 19981 l. Extensive introductions to CLP 
can be found in (IJaffar and Maher 19941 [Mar riott and Stuckey 1998 ). The key notions of 
CLP are those of an algebra and an associated constraint solver over a class of constraints, 
namely a set of first order formulas including the always satisfiable constraint true, the un- 
satisfiable constraint false, and closed under variable renaming, conjunction and existential 
quantification. If c is a constraint, we write 3c for its existential closure. We consider ideal 
CLP(C), i.e., we require the existence of a constraint solver solve mapping in finite time 
each constraint to true or false such that if soJvc(c) = false then the constraint 3c is 
false with respect to C and if solvc(c) = true then the constraint 3c is true with respect to 
C. The associated domain is denoted Dp. Given a constraint c, a solution of c is a mapping 
8 from the set of variables to Dc such that cQ is true with respect to C. The set of predicate 
symbols associated with C is denoted lie. We are interested in the following domains and 
languages: 

• N. The predicate symbols are = and >, the function symbols are 0, 1, and +. 

• Q and K. The predicate and function symbols are as above. Q + and R + restrict Q 
and M to non-negative numbers. 

Given a CLP(C)-program P, we define Yip as the set of user-defined predicate symbols 
appearing in P. We restrict our attention to flat programs, i.e., finite sets of rules in aflat 
form. So each rule is of form: either qoifo) <— c or qo(yo) «— . . . ,q n {y n ) where c 

is a constraint, go, . . . ,q n e Yip, yo,... ,y n denote tuples of distinct variables, f}" = oy~i = ' 
and the set of free variables of the constraint c is included in U" =0 yi. Flat queries are de- 
fined accordingly. A binary program is a flat program such that all rules have no more 
than one user-defined body subgoal. The C-base Bp is defined as {p(di, . . . ,d„) \ p G 
Yip, (di,... ,d n ) € (Dc)"}. For a flat query Q of the form c,A\ , . . . ,A n , the set of ground 
instances of Q, denoted ground c (Q), is the set of conjunctions of the form Ai9, . . . ,A„0 
where is a solution of c. The notion of groundedness is extended to flat rules and pro- 
grams. 

Example 1 

Consider the following CLP(Q) program P: 

r\ p{x) <— x = 2. 

n P (x) <- o = i. 

r 3 p(x) «- 72>x,y =x+l,p(y). 

This program is a binary program, groundq(r\) is {p(2)}, groundq(r2) is 0, ground^ (n) 
is an infinite set that contains, among others, p(12) <— p(73) and p(\/2) <— p(3/2), and 
groundq(P) = groundq(ri) U groundqfa) U groundq^). Note that ground instances do 
not contain any constraint. 

We now discuss the operational semantics of CLP-programs we consider in this paper. 
A state of computation is a pair (A\ , . . . ,A„\\c). We further assume that one of the atoms in 
Ai,... ,A„, say A,-, is selected for resolution by a selection rule. The operational semantics 
can be expressed by means of the following rewriting rules: 
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• (A\, . . , ,A„\\c) rewrites to (n|| false) if there exists a fresh rule Aj <— c',B\ , . . . ,B m in 
P such that c A (A, = A ■) A c' is unsatisfiable; 

• (Ai,...,A n \\c) rewrites to (A h . . . ,Ai-\,B\, . . . ,B m ,A i+h . . . ,A n \\c AA; =A\ Ac') if 
there exists a fresh rule Aj «— c',Bi, . . . ,B m inP such that c A (A,- = Aj) Ac' is satisfi- 
able. 

A derivation from a state So is a finite or infinite sequence of states So, Si,... ,S„, . . . such 
that each Si can be rewritten as 5,-+ 1 . A ground state is a state (A i , . . . , A„ \ true) where each 
A, belongs to Bp. We say that a CLP(C) program P is terminating if every derivation start- 
ing from any ground state via any selection rule is finite, under the operational semantics 
defined above. 

To characterize this notion of termination, we use the notion of level mapping. A level 
mapping for a constraint domain C is a function | • | : Bp — ► R. We adapt the idea of recur- 
rence, originally introduced in (IBezem 19931 1, to CLP: 

Definition 1 

Let P be a fiat CLP(C) program, and • | : C-base — > R be a level mapping. P is called 
recurrent with respect to | ■ | if there exists a real number e > such that, for every A <— 
B U ...,B„ e ground c (P), \A\ e R+, and \Bj\ G R+, |A| > +e for all i, 1 < z < n. We 
say that P is recurrent if there exists a level-mapping such that P is recurrent with respect 
to it. 

Observe that rules of the form p(x) <— c are not taken into account by the definition 
above. Moreover, without loss of generality, we may fix e to 1: if P is recurrent in this 
narrow sense, P is trivially recurrent with respect to Definition Q] Conversely, since e > 0, 
we can safely multiply the values of the level mapping by 1 /e. 

Theorem 1 

(IBezem 19931 1 P is recurrent if and only if P is terminating. 

3 Aim-recurrent programs 

Let us consider programs that can be analyzed by means of affine level mappings. 
Definition 2 

A level mapping | • | is called affine if for any n-ary predicate symbol p £ Yip, there exist 
real numbers p. p j, < i < n, such that for any atom p(e\,...,e„) G Bp 1 : 

n 

\p{ei,...,e„) \ =p l ,. Q + Y J L>p,i e i 

i'=i 

So for a given atom p(e), its affine level mapping is a linear combination of e shifted by 
a constant. We can define the class of programs we are interested in: 

Definition 3 

Let P be a fiat CLP(C) program. We say that P is aim-recurrent if there exists an affine 
level mapping | • | such that P is recurrent with respect to it. 
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Example 2 

The CLP(Q) program P from Example[T]is aim-recurrent with respect to \p(x) \ =73 —x. 

Clearly, if P is aim-recurrent, then P is recurrent thus terminating. Let us show that aim- 
recurrence can be efficiently decided. We start with proving this result for binary programs. 

Theorem 2 

Aim-recurrence of a binary constraint logic program P over Q, Q + , K and M + is decidable 
in polynomial time with respect to the size of P. 

Proof 

The proof is constructive: we provide a decision procedure for aim-recurrence of binary 
constraint logic programs over Q,Q + ,R and R + . The decision procedure extends the al- 
gorithm proposed in ( So hn and Van Gelder 19911 1 for termination of Prolog programs (ab- 
stracted as CLP(N) programs) to binary CLP(C) where C is Q,Q+,M or M+. The algo- 
rithm tries to find an affine level mapping showing that P is aim-recurrent by examining 
each user-defined predicate symbol p of a binary CLP program P in turn (the precise order 
does not matter). For every rule r, say p(x p ) <— c,q(x q ), we test the satisfiability of c. For 
the domains we consider, it can be done in polynomial time ( |Khachiyan 1979 ). If c is not 
satisfiable, we disregard this rule. Otherwise, let n p and n q be the arities of p and q. For the 
rule r, recurrence is equivalent to: 

Ch^ [\p(x P )\ > 1 + \q(x q ) | A \q(x~ q ) | > 0] (1) 

Note that the condition c — » \p(x p ) \ > can be omitted as it is implied by (Q~|). Formula (HJ 
is logically equivalent to C |= c — > \p(x p ) \ > 1 + \q(x q ) | and C |= c — > \q(x q ) \ > 0. Let x p be 
(jc p ,i,. . . ,x Ptnp ), x q be (x 9 j, . . . ,x qtnq ) and let p pfi , . . . ,fi Pt „ p ,p qt o, . . . ,fi q , nq elbe such that 
for any atom p[e\ , . . . , e„ p ) € Bp and any atom q(e u ..., e„ q ) € Bpi \p(e\ , . . . ,e np ) \ = p p _ + 
Y^^p-pjei and \q{ei,. . . ,e nq )\ = p q . + Y^UVqid- Hence, c should imply (p pfi - p qfi ) + 
lH^P-pjXpj + YUli (-p- q ,i)x q j > 1 and fi qfi + I,"liLtqjXq,i > 0. For the sake of uniformity, 
we rewrite the second inequality as p.q,o + Hj=iOxp,i + Y.ililJq,i x q.i ^ 0. Both inequalities 
can be presented using the scalar product notation as p\x > 1 and p'x > 0, where: 

X = (XQiXp^l) . . . iXp^ np ^X qi \ } . . . )X q jiq) 

xq is a new variable fixed to 1 and used to obtain the free coefficient in the product 

ft = (Hp,0—flq,0,Hp,l,---,Hp,n p ,—flq,l,---,—Hg,n q ) 
fl' = {^q,0,0,... ,0,/i 9j l,...,/i ?j „ ? ). 

Hence, the binary rule r gives rise to the following two pseudo linear programming 
problems. The problems are pseudo linear rather than linear because symbolic parameters 
appear in the objective functions. 

minimise 9 = px subject to c Axq = 1 (2) 

minimise 8 = p'x subject to c Axq = 1 (3) 
We note that c Axq = 1 is satisfiable as c is satisfiable and xq is a new variable, and we 



rewrite c Axq = 1 as Ax > b in the standard way ( Schrijver 1986). An affine level mapping 
| ■ | ensuring recurrence exists at least for this rule if and only if 0* > 1 and 8* > 0, where 
0* and 8* denote the minima of the corresponding objective functions. Because of the 
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symbolic constants p p j and fj g j, neither (0 nor ([3]) is a linear programming problem. Now, 
the idea is to consider the dual form: 



maximise r| = b T y subject to A T y = /j T Ay>0 



(4) 



maximise y=bz subject to A z = ji Az > (5) 
where y and z are tuples of adequate length of new variables. By the duality theorem of 
linear programming which holds in C (see ( Schrijver 1986) for instance), we have 8* = T|* 
and 8* = y*. Furthermore, we observe that ft appears linearly in the dual problem (0. 
Hence the constraints of (0 can be rewritten, by adding r| > 1 as a set of linear inequations 
denoted Sr~ l+il . Similarly, the constraints of can be rewritten, by adding y > as a 
set of linear inequations, denoted Sf~°. Let us define defn/>(/?) as the set of binary rules 
defining p in P, S p as the conjunction AredeftWp) [$ - l+q AS?~ ], and Sp as the conjunction 
f\pen P Sp' We have by construction Sp is satisfiable if and only if there exists a affine level 
mapping ensuring recurrence of P. 

Moreover, as P is a finite set of binary rules, computing Sp can be done in polynomial 
time with respect to the size of P and results in a constraint the size of which is also 
polynomial with respect to the size of P. Finally, testing satisfiability of Sp in Q, Q + , R, 
and R + can be done in polynomial time ( Khachiyan 1979 ). □ 



Example 3 

Applying the algorithm to the example Q] we obtain the following two pseudo linear pro- 
gramming problems corresponding to (0 and (0, respectively: 

minimise 8 = /j p axi — H P ,i x 2 subject to 72 > x\ Ax2 = JCi + 1 Axo = 1 

minimise 8 = /u p fi + H Py \X2 subject to 72 > x\ Ax2 =X\ + 1 A*o = 1 
Rewriting the system of constraints as Ax > b and switching to the dual form, we get the 
system Sp: 



■n =yi - 

Tl > 1, 

y\ -yi 
y*-ys 
yi>o, 

J2>0, 

yi >o, 

J4>0, 

ys >0 



-y2-72*y3+y4-ys, 



0,-y 3 



-y4+y5=Hp,i, 



> u < 



Y = Zl - Z2 - 72 * Z3 + Z4 - Z5 , 

y>o, 

Zl —Z2 =flp,0, 
-Z3-Z4+Z5 = 0, 
Z4 —Z5 =Vp,\i 

zi >0, 

Z2 >0, 
Z3>0, 
Z4>0, 
Z5 >0 



Since Sp is satisfiable, P is aim-recurrent. Note that projecting Sp onto the fi p /s gives 
{/^O + 73 * /^i > < — 1}. Any solution to this last constraint is a level mapping 

ensuring aim-recurrence of P. 

An immediate consequence of the result above is that recurrence with affine level map- 
pings is also P-time decidable for non-binary CLP(M) program with rules which contain 
more than one atom in their bodies. Formally, the following theorem holds. 
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Theorem 3 

Aim-recurrence of a constraint logic program P over Q,Q + ,R and R + is decidable in 
polynomial time with respect to the size of P. 

Proof 

Let P be a constraint logic program. Let P' be the binary constraint logic program such 
that for every rule qo(yo) <— c, qi (y\ ),..., q n (/«) with « > 1 in P, P' contains the following 
rules: 

qo(?o) <~ c,qi(fi). 
qaifo) ^c,q n {y n ). 

and nothing else. From Definition Q] we note that P is recurrent if and only if P' is re- 
current. Moreover, the size of P' is polynomial in the size of P. Hence, by Theorems [2] 
aim-recurrence of P' is P-time decidable. □ 

Although the technique above is not complete for programs over N, it is a sound way to 
prove recurrence of programs over this domain: if a program is recurrent over Q, it is also 
recurrent over N. For binary programs, as we allow negative coefficients in the level map- 
ping, we get a more powerful criterion than the one proposed in ( Soh n and Van Gelder 19911 l. 
For instance, termination of Example Q] (considered as a CLP(N) program) cannot be 
proved by (ISohn and Van Gelder 199 11 1. 

For binary CLP(Q) programs, the decision procedure described above has been pro- 
totyped in SICStus Prolog (ISICS 20051 1 using the Simplex algorithm ( |Dantzig 195T) and 



a Fourier-based projection operator (Holzb aur 19951 1 to ease manual verification. There- 
fore, the complexity of the prototype is not polynomial. The implementation is available at 

|http : //www .univ- reunion . f r/~gcc/ soft/blnterm4q.tgz| 



4 Related Works 

The basic idea of identifying decidable and undecidable subsets of logic programs goes 
back to (IDevienne et al. 1993I I. 

Recently, decidability of classes of imperative programs has been studied in dCousot 20051 
Podelski and Rybalchenko 2004|ITiwari 2004t . Tiwari considers real- valued programs with 
no nested loops and no branching inside a loop dTiwari 2004l > . Such programs correspond 
to one-binary -rule CLP(R). The author provides decidability results for subclasses of these 
programs. Our approach does not restrict nesting of loops and it allows internal branching. 
While in general termination of such programs is undecidable (Tiwari 2004), we identified 
a subclass of programs with decidable termination property. Termination of the following 
CLP(M) program and its imperative equivalent can be shown by our method but not by the 
one proposed in dTiwari 2 004). 

Example 4 

q(x) <— — 20 < x,x < 20,y + 5 — x,q(y). 
q(x) <— < x,x < 100, y + 1 = x,q(y). 
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while ((-20 < x < 20) or (0 < x < 100)) do 
if (-20 <x < 20) x = x-5 ti 
if (0<x< 100) x = x-\ fi 

od 

Similarly to (Tiwari 2004), Podelski and Rybalchenko (2004) have considered programs 
with no nested loops and no branching inside a loop. However, they focused on integer 
programs and provide a polynomial time decidability technique for a subclass of such 
programs. In case of general programs their technique can be applied to provide a sufficient 
condition for liveness. 

In a recent paper, Cousot (2005 ) applied abstraction techniques and langrangian relax- 
ation to prove termination. Extension of the basic technique should be able to analyse 
loops with disjunctions in their condition such as Example [4] However, complexity of the 
approach is not discussed and it is not clear whether the technique is complete for some 
class of programs. 

One might like to investigate a more expressive language of constraints including poly- 
nomials. Recall that we require the constraints domain to be ideal, i.e., one needs a decision 
procedure for existentially closed conjunctions. Such a decision procedure exists, for in- 
stance, for real-closed fields such as R (ITarski 19311 |Renegar 1992). For some domains 
such as Q, existence of a decision procedure is still an open problem, although it seems 
to be unlikely (Pheidas 2000). If one restricts attention to real-closed fields, one might 
even consider polynomial level-mappings of a certain power rather than the affine ones. 
One can show that in this case proving recurrence is equivalent to determining satisfia- 
bility of the equivalent quantifier- free formula (Tars ki 19311 ITarski 19511 ). Hence, recur- 
rence is still decidable in this case. Although the known complexity bound of determining 
the equivalent quantifier-free formula given an existential formula is a double exponen- 
tial (IBasu et al. 19961 ICollins 19751 ). to the best of our knowledge the complexity of the 
subclass of formulas which we obtain is an open question. 

5 Conclusion 

In this paper we have considered constraints solving over the rationals and the reals. For 
these domains we have identified a class of CLP programs such that an affine level map- 
ping is sufficient to prove their recurrence. We have seen that membership to this class 
is decidable and presented a polynomial-time decision procedure. The decision procedure 
can also be used as a sound termination proof technique for binary CLP(N) and has been 
prototyped in SICStus Prolog for binary CLP(Q). 
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